Collapse-revival dynamics in strongly laser-driven electrons 
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The relativistic quantum dynamics of an electron in an intense single-mode quantized electro- 
magnetic field is investigated with special emphasis on the spin degree of freedom. In addition to 
fast spin oscillations at the laser frequency, a second time scale is identified due to the intensity 
dependent emissions and absorptions of field quanta. In analogy to the well-known phenomenon in 
atoms at moderate laser intensity, we put forward the conditions of collapses and revivals for the 
spin evolution in laser-driven electrons starting at feasible 10 18 W/cm 2 . 
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The motion of a relativistic electron in the field of a 
plane electromagnetic wave was first evaluated by Volkov 
PQ . While the solution of this problem initiated numer- 
ous investigations of quantum electrodynamical (QED) 
effects in the presence of a strong electromagnetic field, 
the fundamental derivations of this part of QED were de- 
veloped in a series of papers [2H1] about fifty years ago. 
Nevertheless, with the development of powerful lasers 
providing intensities up to I ~ 10 22 W/cm 2 [5] and ambi- 
tions plans beyond [6], the analysis of various quantum 
phenomena [7] such as multi-photon Compton scattering 
[SHIP], electron-positron pair creation [IT], spin effects 
[12] . and quantum plasmas [J3] has recently become in- 
creasingly relevant. In all this research, however, the 
external electromagnetic field is considered as classical, 
and its quantum fluctuations are generally believed to be 
negligible for interactions with very strong fields. 

At the same time, it is well known from quantum op- 
tics that the dynamics of an atom in a comparably weaker 
and resonant laser field depends on the quantum fluctua- 
tions of the field. One of the most interesting phenomena 
of this kind is the collapse-revival effect in the evolution 
of the Jaynes-Cummings model [T31 [15] for a two-level 
atom. This effect was predicted theoretically [16] and 
later observed experimentally [17]. Its qualitative expla- 
nation and analytical description was also given by [T51 - 
[2"U] . It was shown that the evolution of the population 
of the atomic states is characterized by two time scales. 
The first time scale is the period of the Rabi oscillations 
while the second slower one is defined by the collapse and 
revival times of the populations being associated with the 
absorption and emission of the field quanta. 

An electron can on one hand be considered as a two- 
level system with regard to spin space. On the other 
hand the electron has no other internal quantum degrees 
of freedom such as the atom and especially strong laser- 
electron dynamics is generally considered to be perfectly 
described within a classical picture. For that reason, the 
questions of the feasibility of collapse-revival dynamics 
for laser-driven electrons arises. In this Letter, we present 
the electron spin evolution in a single-mode quantized 
field within the framework of Dirac theory. We show 
that an electron may exhibit collapse and revival dynam- 
ics similar to that of two-level atoms and demonstrate a 



possible implementation of this pure quantum effect al- 
ready at widely available laser intensities in the moder- 
ately relativistic regime. 

For an electron in an external electromagnetic field, 
two different time scales are characteristic [21]. One of 
these is defined by the frequency uj of the electromagnetic 
field, i.e. in natural units h — c = 1 as used in the 
complete article: To = 2tt/uj ~ 1.2 • 10~ 4 cm, for the 
typical photon energy oj = leV. The other scale is called 
coherence time T c [5] and is defined by the time required 
to induce a 2-7T phase shift in the Volkov state of the 
electron. 

The existence of these two time scales is expected to 
be simultaneously observable when the total interaction 
time of the electron in the laser field T; satisfies 

T; > T c » 7b, 

which, as shown later, may indeed be fulfilled for realistic 
parameters of the laser pulse. 

In a more general view, the slow oscillations are char- 
acterized by the radiation shift of the electron's quasi- 
energy, which is defined via the Dirac equation with mass 
operator [22] , The shift appears both due to the quan- 
tum nature of the driving electromagnetic field and the 
environmental vacuum field. However, the laser-driven 
processes of emission and absorption of photons with fre- 
quency uj are the main contribution to the transition am- 
plitudes that are proportional to the number of quanta in 
the radiation mode [23] • Therefore, the system of an elec- 
tron and a one mode quantized field can be investigated 
as a closed system in a non-perturbative way. 

We benefit from the fact that the state vector of a 
system with definite total quasi-momentum in a single- 
mode quantized field can be found exactly. For the first 
time, the exact solution of the Dirac equation 

- 7^ - m) V = 0, (1) 

with the potential containing the only significant mode 

m 

A * = ~7rv ( ae ^ (fer) + aV(M ) . ( 2 ) 

Dirac matrices 7^, polarization vector e M , normalization 
volume V, photon four-vector k, photon annihilation and 
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creation operators a and a> of the laser mode, electron 
charge e and mass m, was found by Berson |24j . via the 
Bargmann representation of the creation and annihila- 
tion operators. This solution, however, is not convenient 
for our calculation, and we obtained correspondingly the 
exact solution of Equation ([TJ) in operator form, reading 



1 c -i{gx)+ia^ a(kx) [ kb 
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with / = 7 M / M for any four-vector /, the total four- 
momentum of the system of the electron and the field 
q, b» = ee^/VZVu, a = -(q-b)/((q- k) - 2b 2 ), chrj = 
l/2( A /x+ 1/Vx), x = l/y/l - 2b 2 /{q-k), the quasi- 
momentum of the electron p n , and e n — \Jp 2 l + rn 2 . 
The quasi-momentum of the electron depends on the field 
quantum number, yielding 



p n = q - k\ 
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The set of solutions ^ is orthogonal and complete, i. e. 
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In order to demonstrate that an electron in a single- 
mode quantized field evolves in a similar way as an atomic 
two-level system, we consider the observable value of the 
electron spin: 

(V>|7°7V^O - r'M 



s»{x,t) = 



(6) 



Before the electron crosses the boundary of the light 
beam at t = 0, it is described by the free solution of the 
Dirac equation 



^ = e iPo<r-r „)!« 
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and the field is in the coherent state 
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with the average number of photons n — (3 2 . 

The variable ro defines the position of the electron 
where we have to take into account that in a real exper- 
iment many electrons are distributed in a bunch. As a 
consequence we will later average over this variable by 
integrating over the appropriate density [251 126] . 

During the interaction the system of the electron and 
the field is described by the wave function which is a 
linear combination of the exact solutions (|3j) , i. e. 



*(a;) = J dq^2C q , n ij qn (x); x=(t,x) 



(9) 



\k\ = i 



A = b (ae- l{hx) + a ] V^') 



Figure 1. (Color online) Electron motion in a single-mode 
quantized field without momentum change (coherent scatter- 
ing). p is the initial momentum of the electron, b the polar- 
ization vector of the field, the z axis is directed along the wave 
vector fc, z± denotes the plane perpendicular to z, a and 
are the annihilation and creation operators, respectively, uj is 
the frequency of field, V a normalization volume, and e is the 
electron charge. 



After crossing the border of the field, the interaction 



holds for to 



j, where d is the "thickness" of the 



laser pulse (see Fig. [T]), v the injection velocity of the 
electron, and 9 is the angle between k and Vo. 

At t = to, the interaction is turned off, the electron 
becomes free, and the detector may measure the spin 
s(to). To guarantee that only the spin degree of freedom 
changes, we will analyse the evolution for the polarization 
state of electrons only with the same momentum as in the 
primary beam. This case is essentially different from the 
problem of the polarization of an electron as a result of 
Compton scattering [5]. 

A change in the interaction time to will lead to different 
spin values. This time can be changed in two ways. The 
first is a change in the initial energy of the electron. The 
larger the electron energy is, the less time it spends in 
the field. The second way is a change in the interaction 
time, changing the angle 9. The closer the angle 9 is to 
or 7T, the more time the electron interacts with the field. 

Here, we investigate the case of a strong field such 
that the average number of photons n — f3 2 is large. 
Therefore, the main contributions to the spin value ^ 
will arise from n 3> 1 and the coefficients C qyU can be 
evaluated in this limit. This allows us to calculate the 
average value of the spin analytically. Let us consider 
here the outline of the calculation only. 

From the initial condition 



^ Po (r-r ) U (Po) I 



t=0' 



(10) 



the unknown coefficients C q ^ n in the limit of large n were 
found. Then the wave function ^ was inserted into the 
spin definition ^ and the matrix elements between the 
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field states T n ' <n were calculated. Actually, these matrix 
elements were evaluated using the cumulant method |27| , 
which corresponds to a non-perturbative partial summa- 
tion of the infinite series. The integrals over q, q' , and the 
summation over the field quantum number n were carried 
out analytically. The evaluation of the summation over 
the field quantum number n was performed analogously 
to [IS] by means of the saddle point method. The saddle 
point no has quadratic and linear terms in f3. The zero 
order approximation no = /3 2 gives the quasi-classical 
Volkov limit, and the first order corrections which are 
proportional to /3a induce the quantum fluctuations of 
interest here. In the last step, an averaging over the ini- 
tial electron coordinate ro was carried out and the matrix 
element in spinor space was calculated. As a result, the 
dependence of the spin on time was found in the form 

( S "(t)) = - -fc^K-fc)/ 2 ^ (i + Ren 2 ) 

e e (Po-ky 
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where Ji is the Bessel function of order / and ao an initial 
spin four-vector of the electron. This value of the spin 
may be compared with the time averaged Volkov one, i. 
e. 
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where the bar on top of a symbol denotes the averaging 
in time and A represents the four-potential of the clas- 
sical field. Actually, this averaging corresponds to that 
over the initial phases of the electron in the beam when 
entering the space with the field. 

For the analysis of the evolution of the electron po- 
larization, we also need to calculate the electron current 
density j M (x, t) = (ip\ r y fJ -5(x — r')]^} / (tp\tp) averaged over 
the electron initial coordinate Tq. This value can be eval- 
uated in the same way as considered above for the elec- 
tron spin density, yielding 
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The electron polarization is characterized by its 
helicity — the projection of the spin s onto the direction of 
its momentum v = p a /po- The eigenvalues of the helicity 
operator for the free electron are ±1. In accordance with 
the general rules of quantum mechanics, the expectation 
value of the helicity operator can be written as 



(E-i/)=pi •!+?_! •(-!), 



(14) 



with pi and p_i being the probabilities of observing the 
electron with helicity +1 and —1, respectively. Usually, 



one considers at the beginning the situation of p\ being 
unity and p_i zero, and the interest is how p\ and p-\ 
are modified due to the interaction. 

However, in real experiments, not a single electron, but 
rather a beam of electrons is considered with N being the 
number of electrons in the beam. If the initial electron 
beam was fully polarized, some quantity of electrons with 
opposite polarization should appear after the interaction 
is finished (see Fig. [I]). Namely, the number N_i of such 
electrons would be equal to p-\ N if in the initial moment 
of time t — the helicity of the electron was +1. 

To calculate the final polarization, we choose the co- 
ordinate system with z-axis directed along k, an initial 
electron momentum p displaying an angle 9 with the 
z-axis, and a field being linearly polarized and directed 
perpendicular to the z-axis (Fig. 0. We assume that at 
the initial moment of time the vector £ being the aver- 
age spin in the electron rest frame is directed along the 
momentum in the laboratory frame. As a consequence 
p£ = P, which corresponds to a helicity of the electron 
of+1. 

In order to find the polarization of the electron after 
the interaction with the field, we project the spin vector s 
on the electron momentum direction v = (0,sin#,cos#) 
and normalize it via the probability of finding an electron 
at the observation point, i.e., dividing by j°(t), the zeroth 
component of the current density vector, defined in (13 I. 
Using this fact, Equation (14) and the condition that 
Pi + p—i = 1, the probabilities to find the electron in 
the transmitted and thus non scattered beam in various 
polarization states can be calculated as 



Pi = 



P-i 



111 + e.f - 2£.9 
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(15) 
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where / = («o — cos 9) cos#/(l — v cos6>) 2 (1 + Re n 2 ), 
g = u o sin0/(l - v cos6>)Re IT, f l = (1 + Re II 2 ) /(l - 
^o cos#), and the dimcnsionless parameter £ = /36/eo was 
introduced. 

These probabilities may be compared with the Volkov 
probabilities written in the same variables, i. e. 



1 1 1 + eiv (m , 

P-i = o - o ' (16) 

2 2 



where /„ = (vq — cos#)cos#/(l — Vocos9) 2 and f\ v = 
1/(1 -v cos 6>). 

One can see that contrary to the Volkov case, the fully 
quantum probabilities depend on time. This means that 
the quantum fluctuations can change the system dynam- 
ics. 

Let us find the parameter which governs the slow os- 
cillations in the time evolution in ( 15 ). For this purpose, 



we investigate the system's evolution for small times t. 
In this case, the sine inside the Bessel functions can be 
expanded in its Taylor series and one can find that the 
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there are oscillations which appear due to the quantum 
nature of the electromagnetic field. The average of these 
oscillations represents the quasi-classical Volkov proba- 
bility. 

In real experiments an electron beam involves a spread 
in both initial energy and direction. We ensured that an 
uncertainty of one percent in energy and in the angular 
P-i 
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Figure 2. (Color online) The probability of finding the 
electron with flipped polarization as a function of the an- 
gle 8 for an intensity / = f0 18 W/cm 2 , a pulse duration 0.0005 - 
r = f .5 • 10 _3 cm corresponding to 50 fs, a frequency of the 
photon uj = 7.8 ■ 10 4 cm _1 , an initial probability p_i = 0, and 
7 values of the electron equal to 5 and 10. 




amplitude of the quantum fluctuations is defined via the 
parameter £ which corresponds to x in [3]. 

Modern lasers can reach nowadays high intensities 
[5] up to I0 22 W/cm 2 with a pulse duration of about 
30 fs. For our concrete analysis, we choose as intensity 
I = 10 18 W/cm 2 , as pulse duration r = 1.5 ■ 10~ 3 cm -1 
corresponding to 50 fs, and as photon frequency u> = 
7.8 • 10 4 cm _1 , i. e. a wavelength of 800 nm. 

As was already mentioned, the interaction time can be 
adapted either via the electron's energy or by changing 
the electron's path in the light pulse such as e.g. by vary- 
ing the entrance angle 9. Fig. [2] displays the probability 
of finding the particles with the flipped polarization in 
the electron beam as a function of the equal entrance 
and exit angles 9 between p and k. As can be seen 
from the graphs, when the interaction time increases, cor- 
responding to larger angles, characteristic structures of 
the probabilities appear. In addition to fast oscillations 
at the frequency of field u>, there are slow oscillations 
governed by the parameter £. These oscillations occur 
around the mean value, which corresponds to the quasi- 
classical Volkov case. It should be noted that this special 
structure appears only when the field is considered as a 
quantum object: it can not appear in the quasi-classical 
case. 

Fig. [3] shows the dependence of the probability of find- 
ing an electron with flipped polarization as a function of 
the electron energy for the two values of fly-in angle 9 
equal to 30 and 140 degrees. As in the previous case, 



Figure 3. (Color online) The probability of finding an electron 
with flipped polarization as a function of 7 for an intensity 7 = 
10 18 W/cm 2 , a pulse duration r = 1.5 ■ 10 _3 cm corresponding 
to 50 fs, a frequency of the photon uj — 7.8 ■ 10 cm -1 , an 
initial probability p_i = 0, and values of the fly- in angle 9 
equal to 30 and 140 degrees. 

distribution does not change the displayed probabilities 
in above figures visibly. We emphasise from the order 
of magnitude in the probabilities in Fig. [2] & [3] that 
the number of electrons should be well above 1000 and 
that mutual interactions shall be avoided with appropri- 
ate densities. 

Note that the atomic two-level system possesses an 
analogous behaviour. In this case, the level population is 
considered instead of the electron polarization. The fast 
oscillations in the population inversion essentially depend 
on the Rabi frequency while for an electron in a quantized 
field, it can be linked to the oscillations at the frequency 
of the electromagnetic field uj. The structure of the slow 
oscillations is reminiscent of the "collapse-revival" effect 
for inverted populations, according to which the popula- 
tion inversion of the two-level system vanishes but after 
some time revives again. This effect is purely quantum 
mechanical and can not be found in a quasi-classical anal- 
ysis. We stress in concluding, that the presented quan- 
tum results coincide with those via the quasi-classical 
Volkov solution if the parameter /3a is set equal to zero. 

The authors are grateful to Dr. K. Z. Hatsagortsyan, 
Dr. A. Di Piazza and S. Meurcn for useful discussions. 
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